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^ ■ Abstract 

O 

O ■ The Dual Meissner Effect description of QCD in the confining region 

provides ^ behaviour for the gluon propagator and involves the dual gluon 
mass m as a parameter. This is used in the Schwinger- Dyson equation for the 
quarks in the infrared region to exhibit chiral symmetry breaking for light 

^fn ■ quarks. Using the light quark condensate as input, the dual gluon mass is 

determined and its importance in showing the asymptotic free behaviour of 

K> ! the extrinsic curvature coupling in the rigid QCD string is discussed. 
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Confinement of quarks and gluons and chiral symmetry breaking are 
the two important non-perturbative (infra-red) features of QCD and so any 
model of confinement must necessarily produce chiral symmetry breaking 
nonperturbatively. From the phenemenological quark-antiquark confining 
linear potential, most suggestions of confinement imply -^ behaviour for the 
gluon propagator in the infrared region. Such a linear potential violates the 
cluster property in a local quantum field theory [1]. However, Strocchi [2] 
demonstrated that the proof of the cluster property for QED actually fails 
for QCD upon realizing that the Yang-Mills field strength F" in QCD are not 
observables in contrast to the observability of the electric and magnetic fields 
in QED. One of the authors (R.P) motivated by the suggestion of 'tHooft 
[3] and Mandelstam [4] that monopoles must play a crucial role in the in- 
frared regime of QCD and that of Nambu [5] that dual Meissner effect as a 
mechanism of confinement, described the infrared regime of SU{3)QCD by 
proposing magnetic symmetry condition to select the confining configurations 
of A" [6]. By exactly solving this condition, appropriate SU{3) gauge fields 
relevant to describe the infrared regime were obtained. A dual QCD action 
had been derived and the quantum one-loop corrections generated mass for 
the dual gluons by 'gauge mixing mechanism' [7]. The effective action then 
exhibited condensation of monopoles as in the London theory of Meissner ef- 
fect and the dual monopole corresponded to the Abrikosov flux. Subsequently 
the author (R.P) derived a string representation of the above dual action [8] 
and this was that of a rigid string involving besides the Nambu-Goto term, 
an extrinsic curvature dependent term. A string with extrinsic curvature has 
been proposed earlier by Polyakov [9] and Kleinert [10] to describe QCD and 
has been studied in detail using generalized Gauss Map [11]. The extrinsic 
curvature term has gained importance in view of the result of Kleinert [12] 
who showed that it is generated dynamically by the fluctuations of the string 
world-sheet. Subsequently, Kleinert [13] presented a continuum fleld theory 
of compact f/(l) gauge theory which admits monopoles. In contrast to the 
dual Higgs fleld description of quark conflnement, he [13] used the original 
formulation in terms of the compact U{1) gauge fleld and dualized the ac- 
tion to describe dual Meissner effect. The dual gauge flelds acquire mass 
and this gives rise to linear conflning potential. Later, Polyakov considered 
the compact U{1) gauge theory and obtained rigid string action in the large 
area limit of the Wilson loop. So it is natural that the dual Meissner effect 
description [6] gives the rigid string action [8]. However in [6] and [8], the 



dual gluon mass was not determined. 

It is the purpose of this paper first to demonstrate chiral symmetry break- 
ing in QCD using the confining configurations of [6] that gave the dual Meiss- 
ner effect and then use these results to determine the dual gluon mass m. 
The dual gluon mass is an important parameter in the dual superconductor 
picture of confining region of QCD and its determination will be very useful 
in its own right besides giving the relative strength of the extrinsic curvature 
coupling in the QCD string description. The massive dual gluon propagator 
(Eqn.15 of Ref.6) in Euclidean version has been found to be ^ '"" „ and the 

gluon propagator, using the analysis of Nair and Rosenzweig [15], is —^S^u 
in Feynman gauge. As the above dual gluon propagator has been obtained 
in the confining region of QCD, this gluon propagator is valid only in the in- 
frared region. We will examine the chiral symmetry breaking using the above 
propagator for the gluons and Schwinger-Dyson equation for the quark prop- 
agator in the rainbow approximation. The input to determine the dual gluon 
mass m will be the light quark condensate. 

Schwinger-Dyson equation provides a nonperturbative analysis of the 
quark propagator. A complete solution of this equation is extremely dif- 
ficult and various approximations have been used. For a review, see Roberts 
and Schmidt [16]. The Schwinger-Dyson equation for quark propagator is 

S~\p) = Z,{z^.p + m,) + Z,j^^g'^^^S{q)V^'D^,{p-q), (1) 



where S{p) is the full quark propagator, D^^, is the dressed gluon propagator, 
F'^^ is the dressed quark-gluon vertex, mi, is the bare mass of the quark (which 
will be set equal to zero hereafter), Zi is the quark-gluon vertex renormal- 
ization and Z2 is the quark wave function renormalization. In applying (1) 
to QCD in the infrared region, we invoke the following: (1) In the infrared 
region, the QCD coupling g"^ will be taken to be a constant, following Gri- 
bov [17]. This kind of 'freezing' of g"^ has been pointed out earlier [18,19] 
and discussed by Aguilar, Mihara and Natale [20]. Recently, Alkofer and 
Fischer [21] observed that the nonperturbative strong running coupling re- 
sulting from gluon and ghost propagators possesses an infrared fixed point. 
Although the origin of the infrared fixed point is not known yet, we will 



treat g^ constant in the infrared region of QCD. (2) We will take mt = so 
that the p^-dependent mass term in S{p) at p^ = 0, if not zero, will ensure 
chiral symmetry breaking. (3) We will use rainbow approximation and take 
Y^b = Y^ a.ndZi = Z2 = l. Then writing D^^^ij) - q) = S^'^D^^ij) - g), (1) 
becomes 

S-\p) = ^^.p + g^J^^YS{p-q)YD,,{q). (2) 

The standard procedure is to assume a solution of (2) in terms of functions 
A(p^) and M(p^) so as to have 

where 



pM2(p2)+M2(p2)' 

p'C'{p') + B\p') = {p'A\p') + M\p')}-\ (4) 

Using (3) in (2) and taking the trace over the gamma matrices, we obtain 
the following integral equations for C(p^) and -B(p^). 

B{p') 2 f dS 



9'J^,D;iQ)B{{p-qr), 



p2C2(p2)+52(p2) ^ J (27r 

^ P +9 77rTlC{{p-q) 



p2C2(p2)+52(p2) ^ ' ^ J (27r)4 

{2p^(p -qy-p.(p- q)5nD,u{q). (5) 

The considerations in [6] led to \ behaviour for the gluon propagator in the 
infrared region of QCD and we take D^y{q) in the Feynman gauge as 
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DfM) = --t5,J^'- (6) 



Then the integral equation for B{p^) becomes 



BiP^) _ _^^,^,r AM^^. (7) 



After a change of variable as /c^ = (p — g)^, the right side of (7) becomes 

_4™2„2 r d4fc B{k^) 

and carrying out the angular integration, we have this as 

_o 2 2 2 f k^dk B{k^) 

QTT m g J ^2^^4 (p2_fc2)2 

We shall now consider B{k^) as an analytic function of k^ and use Cauchy- 
Riemann representation 

B(k^) = J-fda-l*^. (8) 

Such a representation is justifiable on the ground that quark fields are not 
observables, the observability criteria being taken as the vanishing of the 
commutator of the BRST charge Qb and the corresponding field [22], the 
physical subspace being annihilated by Qb- As the commutator of Qb with 
quark field is not vanishing, quarks are not observables. This ensures that 
the quark fields are not asymptotic states. Further, the above representation 
has been used recently by Anishetty and Kudtarkar [23] in a theory of QCD 
with string tension. Then, we have 

E(p2) Sn^m'^g^ f k^dk 



p'^C'^ip'^) + B^{p^) 2m J (27r)4 



-oo "(A;2-a)(p2_fc2)2- (9) 



Using Feynman parameterization to write 

1 /•! 2x 

dx- 



{k~p'^Y[k'^ — a) Jo [/c2 — xp"^ — (1 — x)aY' ' 



and carrying out the dk integration, the right side of (9) becomes 



4(27ri)(27r)4 J-oo """^ l"J Jo xp2 + {l-x)a- 



The x-integration is performed to give the right side of (9) as 



where we have used (8) to write the first term. We assume that B{a) is a 
smooth function of a. Then a = p^ is a double pole and a = singularity in 
the logarithm will not contribute. Evaluating the double pole contribution 
at a = p^, we find 

I = a^ (10) 

p2C2(p2)+52(p2) 27r2 ' ^ ' 

and therefore from (4) 

p^a\p^) + mHp^) = ^. (11) 

Thus the unknown functions A(p^) and M{p'^) in the solution of the 
Schwinger-Dyson equation are constrained by (11). This relation obtained in 
the rainbow approximation using ^ propagator for gluons, as follows from 
[6], in the Feynman gauge has not been reported in the literature to the 
best of our knowledge and so constitute a new relation. The validity of this 
relation is only in the infrared region of QCD described by the confining 
configurations of Ref.6. It is interesting however to record that the above 
relation (10,11) also follows from the (5-function behaviour of the confined 
dressed gluon propagator in the non-perturbative region. Such a propagator 
has been suggested by Munczek and Nemirovsky [24] as 

D;'M = 6^\6,^-^-^){2nrg'6\q), (12) 

where ^ is a dimensionful constant. This propagator (12) has been exten- 
sively used by Roberts and his co-workers [25] in their study of Schwinger- 
Dyson equation. Using this propagator in the first equation of (5), we find 

p^C^{p')+B^{p^) = —— = a constant, (13) 



which is same as (10) and so (11) with m^ = Gtc'^Q'^. This comparison, of two 
different gluon propagators giving the same relation, shows the importance 
of (10,11) in the confining region of QCD. We wish to state that while the 
propagator (6) has a field theoretic derivation in [6] and gives rise to linear 
potential, the propagator (12) seems to be an ansatz and will not give rise 
to linear potential. 

In a similar manner, we find C{p^) in (5) becomes a constant and this 
is also the result of using (13). The constancy of C(p^) is valid only in the 
infrared region. This in turn gives A(p^) is | a dimensionless constant. Thus 
(11) becomes 

f.MHp^) ^ ^, (14) 

in the confining low momentum region of QCD. From (14) it follows that 
M(0) = -^ and this shows that chiral symmetry is broken in the confining 
region of QCD described by the model of Ref.6. The result consists of two 
parameters, the dual gluon mass m and the QCD coupling g in the infrared 

2 2 

region. We shall denote ^^^ by D and estimate its value by evaluating the 
quark condensate. The relation between the light quark condensate < qq > 
and M{p'^) for three colors [26,27] is 

^ - ^ - -JL [ rl 2 P^M(p^) 



Using (14), we have 



12 
81DVD 



807r2 



(16) 



Using the lattice QCD estimate [28] of < qq >= -{250MeVf for light 
quarks, we find y/D = 534MeV", which is also M(0). The numerical value 
of M(0) is 0.534GeV which compares well with lattice estimate of ~ O.GGeV 

[29]. 



Now we use the above results to predict a value for the dual gluon mass 

2 2 

m. Since our determination of D using light quark condensate fixes ^|^ 

7 



as 0.28G'eV"^, we need to know g"^ for the determination of m. We use the 
expression of Alkofer and Fischer [21] for a{p'^) which is g'^ /{An) after nor- 
mahzing the parameters to q;(M|) = 0.112. Clearly then m depends upon 
p. In Table. 1 we give the numerical values of m. 



Table. 1 



Dual Gluon Mass Vs a{p'^ 



p(GeV) 


a{p') 


Dual Gluon Mass (GeV) 





2.972 


0.384 


0.1 


2.770 


0.398 


0.2 


2.317 


0.436 


0.3 


1.852 


0.487 


0.4 


1.487 


0.544 


0.5 


1.221 


0.600 


0.6 


1.030 


0.654 


0.7 


0.889 


0.703 


0.8 


0.784 


0.748 


0.9 


0.703 


0.790 


1.0 


0.641 


0.828 



From the Table, we find for g"^ ~ 8, the dual gluon mass is ~ 828MeK, which 
agrees well with the estimate of Baker, Ball and Zachariasen [30] who used 
the above value of g"^. 

To summarize, we have used the \ behaviour of the infrared gluon prop- 
agator obtained from the dual Meissner effect description of QCD in the con- 
fining region [6] in the Schwinger- Dyson equation for the quark propagator in 
the rainbow approximation in Feynman gauge. The chiral symmetry break- 
ing is seen through M(0) 7^ 0. Using the light quark condensate value, the 
numerical value of the combination m^g'^ is determined as ^^^ = 0.28Gel^^. 
We have used the expression of Alkofer and Fischer [21] to evaluate g"^ and 
this depends upon the momentum p. The correspponding dual gluon mass 
values are given in Table. 1. For g"^ = 8 (the value used in Ref.30), the 



dual gluon mass is 828 MeV which agrees with that of Ref.30. One of the 
importances of the dual gluon mass is to observe that in the rigid string 
action for QCD derived in Ref.8 starting from Ref.6, the strength of the ex- 
trinsic curvature term is jg^;^ which becomes in view of the above results 

367r3(a28Gey2) = qtt^ ro^28Ge v^ ) ' where A-1 is a measure of the thickness of 
the QCD-string world sheet. Since the strong coupling a is asymptotically 
free, the extrinsic curvature coupling will also be asymptotically free and this 
conclusion agrees with the results of our earlier work on extrinsic curvature 
couphng [11]. 

By using our infrared gluon propagator, the quark-antiquark potential is 

2 2 2 2 

found to be V{r) = ^|;;§-^ and so the string tension is a = ^f;;!-. Using the 

2 2 

value of ^|;^ = 0.28GeV^, as obtained above, the string tension is estimated 
for light quarks as 0.07GeV^ and this agrees with the recent estimate of 
Weda and Tjon [31]. In the action for the rigid string for QCD in [8] derived 
from the dual Meissner effect description, the extrinsic curvature coupling 
in terms of the dual gluon mass and A is given in the previous paragraph. 
Using a = 0.07GeV^, we find for p = 0.4, 0.6, 0.8 and 1.0 GeV, the extrinsic 
couphng (dimensionless) as 10.24, 2.33, 0.977 and 0.401 respectively. These 
estimates show the importance of the extrinsic curvature coupling in the low 
momentum confining region. 
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